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Abstract 

The tilt modulus of a defective Abrikosov vortex lattice pinned by material line defects is com- 
puted using the boson analogy. It tends to infinity at long wavelength, which yields a Bose glass 
state that is robust to the addition of weak point-pinning centers, and which implies a restoring 
force per vortex line for rigid translations about mechanical equilibrium that is independent of 
magnetic field. It also indicates that the Bose glass state breaks into pieces along the direction of 
the correlated pinning centers if the latter have finite length. The critical current is predicted to 
crossover from two dimensional to three dimensional behavior as a function of sample thickness 
along the correlated pinning centers in such case. That crossover notably can occur at a film 
thickness that is much larger than that expected from point pins of comparable strength. The 
above is compared to the dependence on thickness shown by the critical current in certain films of 
high-temperature superconductors currently being developed for wire technology. 
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INTRODUCTION 



It well known that thin films of high-temperature superconductors exhibit larger critical 
currents than their single-crystal counterparts. Thin films of the high-temperature supercon- 
ductor YBa 2 Cu 3 7 _5 (YBCO) grown by pulsed laser deposition (PLD), which are actively 
being developed for wire technology, achieve critical currents that are a significant fraction 
of the maximum depairing current, for example Evidence exists that lines of dislocations 
that run parallel to the crystalline c axis in PLD- YBCO act as correlated pinning centers for 
vortex lines inside of the superconducting state [2J, and thereby give rise to such high critical 
currents. This is confirmed by the peak observed in the critical current of PLD- YBCO at 
orientations of the c axis aligned parallel to an applied magnetic field, as well as by the dra- 
matic enhancement of the former peak after more material defects in the form of nano-rods 
aHg.edpa.aneHotheca.isa.eaddedafl. 

The microstructure described above for PLD-YBCO films immediately suggests that 
the vortex lattice that emerges from the superconducting state in applied magnetic field 
aligned parallel to the c axis is some form of Bose glass characterized by a divergent tilt 
modulus^. In the limit of a rigid vortex lattice, to be expected at large magnetic fields, 
two-dimensional (2D) collective pinning of vortex lines by the material line defects then 
determines the critical current Q] [3] Q] [9] - Recent theoretical calculations that follow this 
line of reasoning find moderate quantitative agreement with the critical current measured 
in films of PLD-YBCO in c-axis magnetic fields of a few to many kG, at liquid nitrogen 
temperature [10j. A potential problem with the Bose glass hypothesis, however, is that the 
material line defects in PLD-YBCO films can meander, or they can be of relatively short 
length |3|fl. That question is addressed in this paper, where we find that the Bose glass 
breaks up into pieces along the direction of the correlated pinning centers when the effective 
length of the latter is less than the film thickness. In particular, the profile of the critical 
current versus film thickness is predicted to reflect a crossover from two-dimensional to three- 
dimensional (3D) collective pinning of the vortex lines U|. This cross-over can occur at a 
length scale that is notably much larger than that expected from point pins of comparable 
strength [l^ . We also find that the unbroken Bose glass state is robust to the addition of 
weak point pins. Good agreement is achieved between the dependence on thickness shown 
by the critical current in certain films of PLD-YBCO at applied magnetic field [3] and that 



2 



predicted for the true Bose glass state [10|. Last, the effective restoring force per vortex line 
due to a rigid translation of the Bose glass about mechanical equilibrium, which is gauged 
by the Labnscn p_Q, is found to depend o„ ly woaMy on applie d magnetic field. 
This prediction agrees with recent measurements of the microwave surface resistance on 
PLD-YBCO films with nano-rod inclusions [jjj]. 

TILT MODULUS OF BOSE GLASS 

Material line defects in thin enough films of PLD-YBCO can be considered to be per- 
fectly parallel to the c-axis. They notably arrange themselves in a manner that resembles a 
snapshot of a 2D liquid 2], as opposed to a gas[15|[l6|. In particular, such correlated pin- 
ning centers do not show clusters or voids. A defective vortex lattice that assumes a hexatic 
Bose glass state will then occur for external magnetic fields aligned in parallel to such a 
microstructure, in the limit of weak correlated pinning. It is characterized by parallel lines 
of edge dislocation defects that are injected into the pristine vortex lattice in order to relieve 
shear stress due to the correlated pins. The former do not show any intrinsic tendency to 
arrange themselves into grain boundaries, however, due to the absence of clusters and voids 



in the "liquid" arrangement of linear pinning centers. (Cf. refs. [15( and [16|.) This results 
in a vortex lattice whose translational order is destroyed at long range by the isolated lines 
of edge dislocations, but which retains long-range orientational order. Collective pinning 
of the dislocation defects by the correlated pins then results both in an elastic response 
to shear and in a net superfluid density. Theoretical calculations [17| and numerical Monte 
Carlo simulations (see fig. [T]and caption) of the corresponding 2D Coulomb gas ensemble fis| 
confirm this picture. 

The correlation length L c (\) for order along the magnetic field direction is infinite for 
a Bose glass state. In the limit of weak correlated pins, the density of dislocation defects 
that thread the correspondingvortex lattice can then be obtained by applying the theory 



of 2D collective pinning [6j [7 ][8|. Each line of unbound edge dislocations is in one-to-one 
correspondence with a well-ordered bundle of vortex lattice, or Larkin domain, that has 
dimensions -R c (|) x Rc{\) in the directions transverse to the magnetic field. The injection 
of the dislocation lines into the pristine vortex lattice then results in plastic creep of each 



Larkin domain by a Burgers vector 19j . The transverse Larkin scale is hence obtained by 
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minimizing the sum of the elastic energy cost due to the edge dislocations with the gain in 
pinning energy due to the translation of a Larkin domain by an elementary Burgers vector 
of the triangular vortex lattice, b = a&. This yields the estimate^] 

R c {\)- 2 = Cln p {fJc m b)\ (1) 

for the density of Larkin domains, which coincides with the density of lines of unbound 
dislocations. Here n p denotes the density of pinned vortex lines, / p denotes the maximum 
pinning force along a material line defect per unit length, and cqq denotes the elastic shear 
modulus of the pristine vortex lattice. The prefactor above is of order jgj C = n/\n(R c /a' di ) 2 , 
where a' di is the core diameter of a dislocation in the vortex lattice. Consider now the 
limit of weak pinning centers that do not crowd together: f p —>■ and 7tr p ■ <C 1, 
respectively, where denotes the density of material line defects, and where r p denotes 
their range. Simple considerations of probability then yield the identity n p jn^ = ne • vrr p 
between the fraction of occupied pinning centers and the product of the density of vortex 
lines, tib, with the effective area of each pinning center. Substituting it plus the estimate 



c 66 = ($o/8 7I "^l) 2 ^b for the shear modulus of the pristine vortex lattice[21| into Eq. ([I]) 
then yields the result R c (\)~ 2 = (v / 37r/2)C , Q(4/ p r p /e o ) 2 n0 for the density of Larkin domains, 
which depends only weakly on magnetic field. Here, denotes the London penetration 
depth and eq = (Qq/AhXl) 2 is the maximum tension of a flux line in the superconductor. All 
of the above is valid in the 2D collective pinning regime that exists at perpendicular magnetic 
fields beyond the threshold B cp = Cq (v^3/2)(4/ p /eo) 2 ^ > o, in which case many vortex lines are 
tinned by material line defects within a Larkin domain of transverse dimensions R c (\)xR c (\) 



We will now exploit the boson analogy for vortex matter in order to compute the uniform 



tilt modulus of the hexatic Bose glass in the absence of point pinning centers 
to a London model set by the free-energy density 

2 

2 W \ dz 



0. 



It amounts 



SB({r}, z) = £ \ex ( ^ ) + E ^o(r t , r,) + £ V p (r t ) (2) 

i>j i 



for vortex lines located at transverse positions {rj(z)}, at a coordinate z along the field 
direction. Here E\ denotes the tension of an isolated vortex line, while the pair potential 
Vo(r, r') describes the interaction between vortex lines at the same longitudinal coordinate 
z. The energy landscape for the correlated pinning centers is described by the potential 



energy V p (r), which is independent of the coordinate z along the field direction. The ther- 
modynamics of this system in the presence of an external tilt stress u-qsl is then set by the 
partition function 



(n, J P[r,(^)])exp[-(A ;B T)- 1 jf * dz[g B ({r}, z) - J d 2 r j B ■ a]], (3) 



which under periodic boundary conditions, Ti(z+L z ) = r^z), is equivalent to a system of 2D 
bosons. Above, jB(r, z) = [ r ~ r i( z )](dri/dz) is the current density within the boson 

analogy. Observe now that the kernel U^ u (lj) of the uniform electromagnetic response for 
alternating current (AC) is connected to this partition function through the proportionality 
relationship 

Z B [a] ocexp[(A;BT)-V^ia-n-a] (4) 

iuj n 

in the limit that the corresponding uniform tilt stress vanishes, a — > 0. Above, the Matsubara 
frequencies iu n are given by the allowed wavenumbers along the magnetic field, q z = 2nn/L z , 
and V = L x L y L z is the volume of the system. The fact that the tilt stress is given by ubSl 
then yields the identity 

Cw(q z ) = n B /U±(u; n ) (5) 

between the uniform tilt modulus and the uniform AC electromagnetic response of the 
2D Bose glass. The subscript "_L" above is a tag for the pure shear component of the 
electromagnetic kernel (Cf. ref. js]). 

The hexatic Bose glass state is clearly a 2D dielectric insulator within the boson analogy. 
Its electromagnetic response can therefore be modeled by the kernel 



IT ± H = (nn/e^/iu 2 - lu 2 ), (6) 

which is dielectric in the low-frequency limit, and which conserves charge by satisfying the 
oscillator f-sum rule. Above, ujq is the natural frequency of the electric dipole degrees of free- 
dom in the Bose glass. The latter correspond to the 2D Larkin domains in reality, since they 
represent the smallest units of well-ordered vortex lattice that can respond independently 
to an applied force. We then have that the above natural frequency is of order the resonant 
frequency for transverse sound inside a Larkin domain of the 2D Bose glass: uj ~ j/R c (\), 
where 7 = (cqq/cu) 1 ^ 2 is the effective mass anisotropy parameter equal to the transverse 
sound speed within the boson analogy. Here C44 = Ub£i is the tilt modulus due to isolated 
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flux lines. After substitution into Eqs. ([6]) for the AC response, the identity (j3j) then yields 
a divergent tilt modulus for the Bose glass at long wavelength 

Cu{q z ) = n B ii[l + {q z L,)-% (7) 

with a longitudinal scale = Uq 1 that is related to the transverse Larkin length by the 
anisotropic scale transformation ~ -R c ( | )/T- 

Expression (171) for the uniform tilt modulus of a Bose glass is the central result of the 

n 

paper. We shall first extract the Labusch parameter [13] from the singular behavior that it 
shows at long wavelength. In particular, observe that the elastic energy density for a periodic 
tilt of the Bose glass by a displacement u at long wavelength, \C^{q z )qluQ^ acquires a 
contribution of the form |&o M o fr° m this divergence, with k = c 66 /-R c (|) 2 . The latter is 
simply the spring constant per unit volume of the restoring force for a rigid translation of 
the Bose glass state about mechanical equilibrium. Using the estimate c 66 = \e riB for the 



shear modulus of the vortex lattice (2l| yields an effective spring constant per vortex line due 
to 2D collective pinning limited by plastic creep, k /n B , that is given by k p = ^e /R c (\) 2 . 
Notice that k p depends only weakly on magnetic field. The Labusch parameter extracted 
from the microwave surface resistance on PLD-YBCO films with nano-rod inclusions also 
exhibits only a weak dependence on external magnetic field aligned parallel to the nano-rods 



or c-axis, 



to point pins 



n he above should be compared to the corresponding Labusch parameter due 
131 ]. which depends strongly on magnetic field. Indeed, given the conjecture 
k P = (c m /R 2 + c 44 / ' L 2 C ) /n-Q for the Labusch parameter due to 3D collective-pinning implies 
that it decays with increasing magnetic field instead as 1/B 2 in such case. Here we have 
assumed anisotropic scaling, and we have usedpj R c cx B. 

We shall next use Eq. flTJ) for the uniform tilt modulus to test how robust the hexatic Bose 
glass is to the addition of point pinning centers. The hexatic Bose glass shows long-range 
orientational order in all directions (see fig. [1]). The addition of point pins will therefore 
break it up into Larkin domains, of dimensions R' c x R' c x L c transverse and parallel to the 
magnetic induction, that tilt in the transverse direction by a distance of order the size of 



the vortex core, £. Such a break-up then has an e 
pinning energy gain per unit volume that sum tojij] 



astic energy cost per unit volume and a 



Here denotes the shear modulus of the hexatic Bose glass, which can be approximated 
by C66 in the limit of weak correlated pinning, while the corresponding tilt modulus C44 is 
given by expression (J7J) evaluated at wavenumber q z = 1/L C . Also, n' denotes the density of 
point pins, while /q denotes the magnitude of their characteristic force. Minimizing 5u with 
respect to the dimensions of the Larkin domains then yields standard results for these 9J: 
L c = L c (-)[1 + (q z L,)- 2 ] and R' c = R c (-)[1 + feL*)- 2 ] 1 / 2 , where L c {-) = 2c M c m e 'MX and 
R c (-) = ^^cI^cIq 2 ^ 2 /n' f^ 2 are respectively the longitudinal Larkin scale and the transverse 
Larkin scale in the absence of correlated pinning centers. The first equation is quadratic in 
terms of the variable L~ l , and it has a formal solution L^ 1 = (2L c (-)) _1 +[(2L c (-))~ 2 — L^ 2 ] 1 / 2 . 
The hexatic Bose glass is therefore robust to the addition of weak point pins. In particular, 
the longitudinal Larkin scale L c remains divergent for L c (-) > L*/2, which is equivalent to 
the inequality 2 3 / 2 R c (-) > R c (\). 

2D-3D CROSSOVER IN CRITICAL CURRENT 

The critical current density of the above hexatic Bose glass, which is robust to the addition 
of point pins, can be obtained by applying 2D collective pinning theory [h|. All vortex lines 
can be considered to be rigid rods. Balancing the Lorentz force against the collective pinning 
force over a Larkin domain (| yields the identity J c B/c = [(n p f 2 + n' p f 2 ) / R 2 } 1/2 for the 
product of the critical current density along the film, J c , with the perpendicular magnetic 
field, B, aligned parallel to the material line defects[l(|. Here n p and n' p denote the density of 
vortex lines pinned by material line defects and the density of interstitial vortex lines pinned 
by material point defects, while f p and /' are the maximum force exerted on the respective 
vortex lines per unit length. Again, it is important to observe that the critical current 
is limited by plastic creep of the vortex lattice due to slip of the quenched-in lines of edge 



dislocations along their respective glide planes 19j. Minimization of the sum of the elastic and 



pinning energy densities then yields a higher density of Larkin domains in the hexatic Bose 
glass with material point defects added by comparison ([T|): R~ 2 = Cl{n p f 2 + n' p f p )/(cQQb) 2 . 
This reflects the injection of extra lines of edge dislocations that relieve shear stress caused 
by point pins. Substitution in turn yields a critical current density, J c = j c + f c , that has a 
component due to correlated pins set by the identity 

Jc B/c = C n p f 2 /c 66 b, (9) 

7 



and that has a component due to point pins set by the ratio j' c /j c = n ' p fp / n pfp- The 
critical current density notably varies as J c oc B- 1 / 2 with magnetic field in the limit of weak 
pinning lOj • 



Consider now a film geometry of thickness r along the axis of the material line defects. 
The forces due to point pins add up statistically along a rigid interstitial vortex line. The 
effective pinning force per unit length experienced by an interstitial vortex line is then given 
f P = fU ( r p r ) at fihn thicknesses r that are much greater than the average separation 
Tp between such pins along the field direction. Again, f' Q denotes the maximum force exerted 
by a point pin. The relative contribution by point pins to the critical current density is then 
predicted to show an inverse dependence on film thickness, j' c /j c = t /t, that is set by the 
scale t = (^p/^ P )(/o 2 //pTp). We thereby obtain the linear dependence on film thickness 
I^ 2D ^ = (to + r)j c for the net critical current per unit width, tJ c . Last, comparison of Eq. 
(JTJ) with Eq. (Q yields the useful expression 

L* = 7-H(3 B/ V2 7/a Ci>)(£ • «vx)(jo/j c )] 1/2 (10) 
for the longitudinal scale characteristic of the Bose glass as a function of the bulk critical 

—1/2 

current density, j c . Here a vx — n B ' is the average distance between vortex lines and 
jo is the depairing current density (see ref. Figure [2] shows a fit to data for the 

critical current versus thickness obtained from a thin film of PLD-YBCO at liquid-nitrogen 



temperature in 1 T magnetic field aligned parallel to the c axis 



121 ] . A bulk critical current 



density j c = 0.22MA/cm 2 is extracted from it. Using a value of jo — 36MA/cm 2 for the 



depairing current (ref. [121]). of £ = llnm for the coherence length, of 7 = 7 for the mass 
anisotropy parameter, and setting Cq = 1 yields a longitudinal scale L* = 24 nm. 

Finally, consider again an arrangement of material line defects that are aligned along the 
c axis, that show no voids or clusters in the transverse directions, but that are broken up 
into relatively long rods of length L ^> L*. It can be realized by meandering material line 
defects [l] (2], where L is the correlation length for alignment along the c-axis, or by artificial 
nano-rod defects^ {4}]. Expression ((Tj) for the divergent tilt modulus indicates proximity to 
the regime of 2D collective pinning, where both L c and L$ are divergent. It hence indicates 
that L c is also large compared to L* here. Second, observe that ~ R c {\)/^ coincides 
with the longitudinal Larkin scale if the rods are considered to be point defects. The 
previous ultimately implies the chain of inequalities L > L c ^> L*. They are consistent 
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with a broken Bose glass state for the vortex lattice, which is threaded by isolated lines 
of edge dislocations of length L c along the c-axis that are connected together by lines of 
screw-dislocations 19] along the transverse directions As in the limiting case of the true 
Bose glass state (fig. [T]), the lines of edge dislocations are injected into the pristine vortex 
lattice in order to relieve shear stress due to the correlated pins. Larkin domains hence 



are finite volumes 



6] 9]. In contrast to their transverse dimensions, however, Larkin domains 



exhibit well defined boundaries along the direction parallel to the correlated pinning centers, 
across which the vortex lattice slips by a Burger's vector due to the presence of the screw 
dislo cations [19 ] . The critical current density expected from this peculiar example of 3D 
collective pinning therefore coincides with that of a Bose glass of thickness L c : J c = j c +j' c , 
with a bulk component due to interstitial vortex lines j' c /j c = t /L c . The critical current 
per unit width then varies with film thickness as /^ 3D - ) = tJ c , showing no offset. Figure [2] 
depicts the predicted dependence on film thickness for the critical current of such a broken 



Bose glass. It notably exhibits 2D-3D cross-over at film thicknesses in the vicinity of L c 11]. 



CONCLUDING REMARKS 

The dependence of the critical current on thickness r shown by certain films of PLD- 
YBCO is in fact consistent with dimensional cross-over at r ~ 1 fim. A previous attempt 
by Gurevich to account for such behavior by collective pinning of individual vortex lines at 



point defects yields a longitudinal Larkin scale L c ~ 10 nm that is too small, however [12]. 
We find here, on the other hand, that pinning due to correlated material defects of length L 
yields a Larkin scale L c that is much longer than that [L* = 24 nm, see Eq. (I10p ] expected 
from point pins of comparable strength if L is much longer than that scale as well. Indeed, 
we predict here that the film thickness at which the critical current crosses over from 2D 
to 3D behavior is of order the effective length of the correlated pinning centers when that 
length satisfies L ^> L*. 

The author thanks Leonardo Civale, Chandan Das Gupta and Sang-il Kim for discussions. 
This work was supported in part by the US Air Force Office of Scientific Research under 
grant no. FA9550-06-1-0479. 
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FIG. 1: Shown is a Delaunay triangulation of a low-temperature groundstate made up of 2016 
vortices that interact logarithmically over a 336 x 336 grid with periodic boundary conditions, 
and that experience an equal number of identically weak 5-function pinning centers arranged in 
a "liquid" fashion. (See ref. [3].) The state shows macroscopic phase coherence and long-range 
hexatic order, with a supefluid density and a hexatic order parameter, respectively, that are 29% 
and 58% of the maximum possible values attained by the perfect triangular vortex lattice. It was 
obtained after simulated annealing from the liquid state down to low temperature, which resulted 
in 371 pinned vortices. A red and green pair of disclinations forms a dislocation (see ref. [ijj]). 
The above Monte Carlo simulation results are consistent with theoretical predictions of a hexatic 
vortex glass state in two dimensions, in the zero-temperature limit, for pinning arrangements that 
do not show any clusters or voids (ref. Josephson coupling between layers then produces a 

3ose glass transition above zero temperature as long as the 2D glass transtion is second order (ref. 
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film thickness (microns) 



FIG. 2: Shown is the dependence on film thickness predicted for the critical current of a defective 
vortex lattice found in a broken Bose glass state (solid line). The dashed and dotted lines are 
extrapolated from the 2D and 3D behaviors, respectively. Measurements made by Sang Kim 
(circles, ref. fl^ ) of the critical current on a thin film of PLD-YBCO at liquid nitrogen temperature 
subject to 1 T magnetic field aligned parallel to the c-axis are fit to the straight dashed line 
predicted by 2D collective pinning (ref. [hJ). This yields an intercept — tq = — 69 nm and a 
slope j c = 0.22MA/cm 2 . Although the value of L c shown here is indeed larger than the lower 
bound = 24 nm [see Eq. ([lUp] and is consistent with the fit to 2D collective pinning, it is only 
hypothetical. 
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